Quantum oscillations in the anomalous phase in Sr3Ru207 
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We report measurements of quantum oscillations detected in the putative nematic phase of Sr3Ru207. Sig- 
nificant improvements in sample purity enabled the resolution of small amplitude dHvA oscillations between 
two first order metamagnetic transitions delimiting the phase. Two distinct frequencies were observed, and their 
amplitudes follow the normal Lifshitz-Kosevich profile. The Fermi surface sheets seem to correspond to a sub- 
set of those detected outside the phase. Variations of the dHvA frequencies are explained in terms of a chemical 
potential shift produced by reaching a peak in the density of states, and an anomalous field dependence of the 
oscillatory amplitude provides information on domains. 



In recent years, the possibiUty that electronic Uquid crys- 
talline states might exist has generated considerable inter- 
est. In such systems, the electron fluid is thought to develop 
textures leading to anisotropics which are not compatible with 
the space group symmetry of the host crystal. Experimen- 
tally, nematic-like transport properties have been observed in 
high-purity two dimensional electron gases l^-Hl, and sig- 
natures of electronic anisotropy have also been reported in 
YBa2Cu307_8 |6, 7]. Combined with observations on the 
subject of this letter, Sr3Ru207, these experiments have stim- 
ulated a body of further theoretical work on the issue ll8l- [l9ll . 

Sr3Ru207 offers an excellent opportunity for experimental 
study of this novel behaviour. In the cleanest single crystals, 
with mean free paths of several thousand angstroms, a well- 
defined phase displaying nematic-like transport anisotropics 
exists below 1 .2 K for magnetic fields between approximately 
7.9 and 8.1 T i^]. The combination of long mean free 
paths and the existence of large single crystals means that it is 
possible to obtain thermodynamic as well as transport infor- 
mation about this phase and its formation ll23ll . 

In this letter, we present a study of the nematic phase 
in Sr3Ru207 using one of the best-established microscopic 
probes of metals, the de Haas - van Alphen (dHvA) effect. 
Performing the project required the growth of a new gener- 
ation of high purity single crystals, screened in a series of 
dHvA measurements to provide the largest oscillatory signals. 
This enabled the observation of oscillations both within the 
nematic phase itself and in its immediate vicinity. 

The crystals were grown using methods published pre- 
viously I24i1 . Measurements with a noise floor of around 
30 pV/VHz were performed in Cambridge using a dilution 
refrigerator. An improvement in signal due to crystal purity 
of a factor of approximately 25 was obtained compared to the 
previous dHv A study 125142711 . 

In Fig. [T|we present plots of the dHvA oscillations, inside 
and outside the nematic phase, with the field oriented along 



the c-axis. Second harmonic detection was used, in order to 
reduce the large non-oscillatory component of the suscepti- 
bility near the metamagnetic transition. Consequently, meta- 
magnetic jumps in the magnetisation appear as sharp asym- 
metric peaks. The top panel of Fig. [T] shows the dHvA sig- 
nal in the region neighbouring the metamagnetic transition, 
where three such peaks can be observed, denoted A, B and 
C. The broad peak A located at 7.5 T corresponds to a meta- 
magnetic cross-over previously observed in several properties 
II20I1 . while the sharp features B and C correspond to metamag- 
netic phase transitions at 7.9 and 8.1 T, bounding the nematic 
phase. 

In the bottom left panel of Fig. [T] we present a blow-up 
of the second harmonic susceptibility signal between 7.9 and 
8.1 T at 15 different temperatures between 58 and 500 mK, 
offset for clarity. Small amplitude oscillations are observed on 
top of a smooth magnetic background, which was modelled 
by fitting a fifth order polynomial (red line) to the data be- 
tween 7.915 and 8.080 T. Polynomials of lower order did not 
reproduce the magnetic background well, while higher orders 
removed a component of the oscillations. In Fig.[Tl right, we 
show the oscillations remaining after the background removal. 
These are suppressed as the temperature increases to 500 mK, 
due to the Lifshitz-Kosevich (LK) amplitude reduction fac- 
tor |28]. The quality of the LK fits suggests that both orbits 
result from Landau quasiparticles. The oscillations were re- 
produced on three samples over four different dHvA runs. 
They were moreover found to be suppressed by rotating the 
field as little as 5° from c-axis. 

The oscillatory data taken in the nematic phase revealed two 
distinct frequencies of 1.0 and 2.3 kT (inset to Fig. 01, with 
masses of 6. 1 ±0.3 nie and 7.4±0.2 nie respectively i29ll . For 
reference, we show, in the top panel of Fig.|2] a spectrum ob- 
tained outside the nematic phase, between 5 and 6.5 T. Peaks 
are observed at 0.43, 0.90, 1.78 and 4.13 kT termed 5, yi, «! 
and a2 respectively Isoll : higher resolution studies of the low 
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FIG. 1: dHvA oscillations in the nematic phase. Top dHvA oscilla- 
tions in d^M/dH^ in the region surrounding the metamagnetic tran- 
sition. From left to right, three magnetic features are visible, denoted 
A, B and C on the diagram (see text). Left The signal in the nematic 
phase region, along with fifth order polynomial background subtrac- 
tion, shown in red. Right Nematic phase oscillations remaining after 
background removal at 15 different temperatures. In the left and right 
graphs, curves are offset for clarity, and include a total amplification 
factor of 10^ from low temperature transformers and preamplifiers. 



FIG. 2: Top Fourier spectrum of the dHvA data on the low field 
side of the transition, between 5 and 6.5 T, with the field aligned 
with the c-axis. Inset Fourier spectra of the nematic phase data at all 
temperatures. Main panel Temperature dependence of the amplitude 
of the nematic phase oscillations. Solid lines represent LK fits to 
the data, which revealed quasiparticle masses of 6.1 ± 0.3 and 
7.4 ±0.2 me for the 1.0 and 2.3 kT peaks respectively. Data plotted 
with crosses correspond to temperatures below 90 mK which were 
corrected, as explained in V29A . and were not used in the non-linear 
fits. 



frequency region reveal further frequencies, termed P and 72, 
at 0.15 kT and 0.11 kT. The associated quasiparticle masses 
rule out the possibility that any of these frequencies is a har- 
monic of any other (see lEsi IztI ISQI for details). The resolu- 
tion with which we can determine the frequency of the 1 .0 kT 
peak in the nematic region is limited by the number of oscilla- 
tions present in the narrow field range, and it seems reasonable 
to associate it with the 0.9 kT peak seen on the low field side. 
At first sight, it is difficult to relate the 2.3 kT peak with those 
seen at low fields, but, as shown in Fig. [3] the ai and 0.2 fre- 
quencies show a strong field dependence in the metamagnetic 
region, and the 2.3 kT frequency is compatible with an inter- 
polation into the nematic phase of the field dependence of the 
ai frequency. This identification is also compatible with mass 
analysis (presented in detail in ref. (27*]) which shows that the 
masses associated with the frequencies in the nematic phase 
are similar to those for the p and ai frequencies (6.6±1 
and 7.0±0.5 nie respectively) outside it. 

The dependence of the frequency F of dHvA oscillations 
in the vicinity of metamagnetism encodes information about 
the mechanism of the transition. As pointed out by Julian 
and co-workers [31], the measured dHvA frequency F{B) oc 
A{B) — BdA{B)/dB, where A{B) is the Fermi surface cross- 
sectional area as a function of magnetic field B. The situation 
most usually considered, sketched in Fig.llJj, is that where the 



coupling between field and magnetisation is field dependent 
(e.g. the effective g factor is field dependent), and leads to 
symmetrical changes in the cross-sectional area of the Fermi 
surface for both spins around the zero field value. This pro- 
duces peaks in the measured frequency of opposite sign. 

The situation in Fig. |4}7 originates from a pronounced 
change of gradient in the dispersion of the band being studied, 
which may appear either below or above the Fermi level. Such 
a dispersion leads to changes in mainly one of the spin split 
cross-sectional areas, and produces a peak in one of the mea- 
sured frequencies, leaving the other almost constant. Large 
frequency splitting is therefore expected through the transi- 
tion region and at fields above it. 

Finally, the third scenario depicted in Fig. ^ is that in 
which a change in chemical potential occurs at a certain 
field value. This would coiTespond to a sudden dispersion 
change in one of the other bands in a multi-band material 
like Sr3Ru207, leading to a redistribution of carriers. This 
chemical potential shift leads to an acceleration of the split- 
ting of the cross-sectional area for one spin, and to a slowdown 
for the other. Consequently, peaks of the same sign for both 
spins are observed in the measured frequencies, and both sat- 
urate at high-field values reduced from those on the low-field 
side [3j]. 
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FIG. 3: Field dependence of the dHvA frequencies situated near 1.8 
and 4.1 kT in the low field side. These feature non-linear changes 
in the metamagnetic region. The data points correspond to the loca- 
tion of peaks in FFTs taken over field windows of equal inverse field 
width of 0.005 T^', and the field value for each data point corre- 
sponds to the inverse of the average inverse field. The red lines rep- 
resent the nematic phase field boundaries, and the red cross the value 
of frequency measured inside the nematic phase. The top panel is 
a blow-up of the data within the green square in the bottom panel. 
Three peaks in F(B), denoted A-C are observed near each of the 
metamagnetic features shown in Fig. [T] Corrections to field values 
due to the intrinsic magnetisation were found to be small. 



FIG. 4: Three scenaria where field dependent dHvA frequencies may 
be encountered. In each case, the dispersion is shown on the left, the 
Fermi surface cross-sectional area as a function of magnetic field, 
A{B), in the centre and the dHvA frequency as a function of magnetic 
field, F{B), on the right. Large changes in the slope of A(6) lead to 
peaks in F{B) [31]. a) A field dependent g factor leads to non-linear 
changes inA(S) for each spin species, producing peaks of opposite 
sign in F(6). b) A flattening of the dispersion leads to a large change 
in F{B) for only one spin species, c) A change in chemical potential 
leads to an acceleration of the change in A{B) for one spin species 
and a deceleration for the other. This leads to peaking of F{B) that 
is similar for both spin species. 



Only the third scenario is compatible with our observations. 
We see peaking of both the ai and a2 frequencies through 
the metamagnetic region, followed by saturation at constant 
values of approximately 1.6 and 3.9 kT respectively, but do 
not resolve any frequency splitting through the transition (see 
Figs. OandgJ;). 

This finding is consistent with the results of Tamai et al. 
ifioll . who showed that a saddle point in the dispersion of Y2 is 
present below the Fermi level (^f)- This feature of the band 
structure leads to a peak in the density of states which is situ- 
ated at around 3 meV below E^. Conservation of the number 
of electrons requires that when approaching this saddle point 
with Zeeman splitting, the hole filling of the peak in the den- 
sity of states by one of the spin species leads to an increase of 
the chemical potential compared to its zero field value. This 
corresponds to an exchange of holes from existing bands to- 
wards the region of the Brillouin zone with a large and ex- 
panding density of states, and produces the observed increase 
in magnetisation. Both ai and 0.2 correspond to hole sheets 
of the Fermi surface, so this is consistent both with the reduc- 
tion of their dHvA frequencies at high fields (above 8. 1 T) and 
with the positive sign of the peaks in the frequencies through 



the metamagnetic region. 

This observation has interesting implications for under- 
standing both the metamagnetism and the effects that meta- 
magnetic fluctuations have on different bands near Ef. The 
(X\ and (X2 sheets are based on the d^z and dy^ orbitals of Ru 
and hence have a more 1-D character than the other Fermi sur- 
face sheets. It appears that these bands contribute only weakly 
to the magnetic moment change, but still couple to the meta- 
magnetism via the charge transfer that changes the chemical 
potential. This in turn implies a sensitivity of the ai and 0.2 
sheets to the fluctuations that occur as one or more other bands 
(most likely 72) approach metamagnetic criticality. 

Having observed well defined orbits in the nematic region 
we can analyze changes in the amplitude of each frequency to 
investigate microscopic scattering on different Fermi surface 
sheets. We first consider the normal region; in the top panel of 
Fig.|2] the relative amplitude of the p to the ai and a.2 signals 
is 1 : 12 : 1.3. The mean free path determined from a Dingle 
analysis that describes damping of quantum oscillations by a 
homogeneous distribution of impurities is 270 nm, fully con- 
sistent with that deduced from the resistivity using the Fermi 
surface from ref. |30]. Going into the nematic region, a naive 
analysis of the resistivity would suggest that the mean free 
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path halves l|20|]. Using the Dingle factor and taking into ac- 
count the Bessel function factor associated with the field mod- 



ulation technique II33I1 we can calculate relative amplitudes for 
a reduced mean free path of 135 nm by extrapolating the am- 
plitudes from just outside this region. This gives an amplitude 
ratio 1:8: 0.5. The observed ratio (inset Fig.|2]main panel) 
is very different: 1 : 0.4 : 0. This strongly indicates that the 
assumption of uniform impurity scattering is no longer valid 
in the nematic region. 

To explain this discrepancy we consider the existence of 
domains as discussed in refs. 112 ill and ll22ll . In those stud- 
ies we commented that the rise in resistivity in the nematic 
region for fields applied parallel to the c-axis and the trans- 
port anisotropy seen in slightly tilted fields were both con- 
sistent with strong scattering at domain walls. For quantum 
oscillations this hypothesis implies that only frequencies due 
to orbits entirely within a single domain survive; therefore, 
the observed orbits can be used as real space callipers of the 
domain's smallest linear dimension. The cyclotron diameters 
of the p, ai and a2 orbits are 270 nm, 385 nm and 580 nm 
respectively ll3411 . so the strong frequency dependence of the 
signal suppression suggests that the sample contains a distri- 
bution of domains with an average wall separation of order 
500 nm. Estimating their linear dimension is an important 
step forward, and should motivate studies with probes capa- 
ble of imaging them directly in real space. 

The above analysis gives clues about the microscopic na- 
ture and mechanism of formation of the novel phase in 
Sr3Ru207, but the most important finding presented here is the 
discovery that quantum oscillations can be observed within it. 
Combined with the observation that the frequencies change 
across the metamagnetic region (Fig.|3]l and recent Hall effect 
measurements showing that the Hall number does not drop 
on entry to the phase [35], our work provides strong evidence 
that the phase diagram of Sr3Ru207 is characterised by a se- 
ries of fluid-fluid transitions. Further, the fluid in the nematic 
phase contains Landau Fermi liquid quasiparticles (Fig. |2j 
whose quantum oscillatory amplitude conforms to the predic- 
tion of the Lifshitz-Kosevich theory. Since dHvA is sensi- 
tive to Fermi surface areas and not to shape, this set of ob- 
servations is compatible with any microscopic model postu- 
lating that the root of the nematic transport properties lies in a 
Pomeranchuk-like Fermi surface distortion. 

In conclusion, we have succeeded in detecting quantum os- 
cillations in the nematic phase of Sr3Ru207. Our results es- 
tablish that at least some of the mobile charges in the phase 
are Landau quasiparticles, and also reveal the importance of a 
chemical potential shift and domain physics in understanding 
the phase's formation and properties. 
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